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ABSTRACT 

Using theoretical formulation of Chung -In Um et al., the ground 
state energy of two-dimensional liquid 4 He interacting through 
Lenard-Jones and Aziz potential were calculated at various 
densities. Our theoretical result indicates that ground state energy 
of two dimensional liquid 4 He interacting with Aziz potential is 
close to GFMC calculation. 

Keywords: Lenard -Jones potential, Aziz potential, GFMC 
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1. INTRODUCTION 

For many years, physicist have used 
hyper netted (HNC) and Fermi hyper netted 
chain (FHNC) theories to investigate the 
properties of highly correlated interacting 
boson and fermion systems such as liquid 
4 He and normal liquid 3 He 12 .The 
HNC/FHNC equations with Jastrow ground- 
state wave functions treat consistently both 
long and short correlations of the system and 
give agreements to some degree with 
experimental results. In spite of qualitative 
successes of the Jastrow ground-state 



function, there exist quantitative differences 
between HNC/FHNC results and 
experimental ones, especially the 
equilibrium ground-state energy, and this 
fact leads one to consider the contributions 
arising from the three-body correlation 
functions. 3 It is well known that the three- 
body correlations do not affect seriously the 
structural properties of ground states, such 
as the radial distribution function and liquid 
structure function, but contribute 
significantly to the ground-state energy. 
Since liquid 3 He system is less dense than 
liquid 4 He, that is, the equilibrium density of 
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the former is 0.0166 A 1 while that of the 
latter is 0.02185 A" 3 in a three-dimensional 
system, effects of three-body correlations in 
3 He are rather small compared to those in 
4 He. At equilibrium density, they make up 
about half the difference between the 
Jastrow results and experimental energy in a 
liquid 3 He system 5 and lower the Jastrow 
ground-state energy by about 10 % in liquid 
3 He. 

A trial ground-state wave function 
including three-body correlation functions in 
addition to two-body correlation effects can 
be constructed in the form 

m . h r H ) = TJ f 3 (r„ ) fl h (»» . r fl . r u )l*> = FIO) 

i<j i<j<k 

Where fifcj) is the two-body 
correlation function which describes spatial 
correlations and depends only upon the 
spatial distance ry between two particle ^(ry, 
fjk, r ki ) is a three-body correlation function, 

and '^is the anti-symmetric product of 
non-interacting single-particle plane waves 
normalized in the system's area Q with N 
fermions, which becomes simply unity for 
boson systems. 

In addition to the three-body 
correlation functions in the ground-state 
wave function, one should take into account 
the contributions coming from the 
elementary diagrams as well as effects of the 
HNC/FHNC equations to obtain better 
results. Unfortunately, there does not exisht 
a general closed expression for the 
elementary diagrams and therefore one must 
calculate each diagram individually, which 
is very time consuming. During the last two 
decades, two approximations for the 
estimation of the elementary diagrams have 



been developed. One is known as the scaling 
approximation 4 and the other is the 
interpolating equation approximation 5 .Both 
techniques have shown their efficiency for 
central two-body correlations, providing 
similar results for the total ground-state 
energy to those of the variational Monte 
Carlo simulation. The scaling approximation 
uses the fact that successive contributions of 
the five, six and higher order elementary 
diagrams are approximately proportional to 
and have very similar spatial behaviour to 
the body elementary diagrams, 4 so that one 
can represent the total elementary diagram 
contribution by means of the readily 
calculable four-body elementary diagrams. 
The interpolating equation approximation 
determines a parameter such that it gives the 
same value for the classical isothermal 
compressibility obtained from both the 
HNC/FHNC equations and Percus-Yevick 6 
equations using the classical pressure 
derivative or the compressibility integral. 

Although large number of 
theoretical works regarding 4 He and 3 He 
system have been reported 7 11 but very few 
of them have been used for two-dimensional 
liquid. 

Chung-In Um et al. have developed 
a theoretical model to study two- 
dimensional quantum fluid. By taking their 
model, one has evaluated the radial 
distribution function g(r) as a function of r 
for two dimensional liquid 4He interacting 
through Aziz potential for various densities. 
In section 2, we have given the mathematical 
formulae used in the evaluation. The 
discussion of the result has been given in 
section 3 of this paper. In the last portion, 
we have shown the tables of the results and 
the reference taken. 
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2. MATHEMATICAL FORMULAE 
USED IN THE EVALUATION 

The Hamiltonian for the considered 
system is given by 



tf = r+v = X7; + ^=-^V:+l4<-ol) 

( (< j "* i i< j 

(1) 



Where V(r) is the interatomic 
interaction. A popular potential theoretical 
investigations is that of Lenard- Jones, which 
has parameters, i.e., the hard-core radius and 



the well depth, and can be written as 



V u (r) = 4e 







12 




6 























,£= \Q22K,a= 2.556 A 
(2) 

However, V u (r) includes only the 
dipole-dipole interaction and does not take 
into account multiple interactions. A more 
physically realistic potential which accounts 
for the self-consistent field Hartree-Fock 
repulsion and multiple interactions is the 
HFDHE 2 potential of Aziz et al. 13 presented 
as 



V u (r) = £ 



Aexp 



-a- 



'm J 



(r 

m 

V r 



(r 

m 

V r 



+ c 



10 



(r 

m 

V r 



10 



F(r) 



(3) 



Where 



F(r) = 



|exp[-(Dr m /r-l) 2 ] 



rlr m <D 
rlr m <D 
(4) 



The values of the constants are 
A = 0.54485 x 10 6 , e*/k B = 10.8 
a = 13.353384, C 6 = 1.3732412 

C 8 = 0.4253785, C 10 = 0.178100 

D= 1.241314, r m = 2.9673 A 



It is well known from the Green's 
function Monte Carlo simulations and other 
variational calculations that the Aziz 
potential gives closer results to experiments 
than V u in three dimensions. One uses both 
potentials and com pare the results with 
those in two dimensions. 

The n-particle distribution function 
is represented as 



P n = (.ii, r n ) = v"N(N-l)x 



JV + (n> OVin, r n )dr n+v ...dr N ) 



(5) 



Where 1 denotes the spatial 
coordinates of the i* particle and v is the 
spin degeneracy of the system (2 for this 
system). In homogeneous system, as N — > go 



and Q — > oo, the single-particle distribution 
function reduces to the density of the system 
so that the radial distribution function can be 
expressed in terms of the density and two- 
particle distribution function P 2 (r) as 
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g(r)=\p 2 (r) = N(N-l) 

P 



yV JVg, rJViA, r n )dr v ...dr N ) 

VP) 



(6) 



One uses the FHNC approximation 
to sum the diagrams arising from the cluster 
properties of g(r). In this scheme, g(r) can be 
decomposed as 



8(r ij ) = g dd (r y ) + 2g de (r & . ) + g ee (r.. ) 



(7) 



Where dd, de and ee represent terms in 
which both i and j are not exchanged, only j 
is exchanged, and both i and j are 
exchanged, respectively. The components of 
g(r;j) are given by 



gdd(rij) = exp[u 2 (ry) + N dd (ry) + E^)] 
gde(rij) = gddfcjX Ndefcj) + E de (rij)] 
gee(rij) = g d d(r i j)[-L 2 (r ij ) / v + N ee (r ij ) 

+ E^y) + (Neefcj) + E ed ( riJ )) 2 ] 

gee (ry) = gdd(rij)L(ry) / v 



(8) 



where 

L(r y ) = -l(k F r y ) + v(N ee (r ij ) + E ee ( rjj )) 
l(x) = 2j 1 (x)/x 



(9) 
(10) 



k F is the Fermi momentum of the system, 
and ji (x) is the Bessel's function of the first 
kind of order 1. N mm represent sums of the 
diagrams and Enm sums of the elementary 
diagrams. The equation g ee (rij) denotes the 
terms in which both i and j are exchanged in 
incomplete exchange loop, and N ee (rjj) and 
Eee(rij) are sums of the nodal diagrams and 
the elementary diagrams in which i and j 
belong to the same permutation loop, 
respectively. 



With the FHNC equation and 
expression for the g(r), one now calculates 
the energy per particle of the system 



1 (w\H\w 



N N N(y/\y/ 



N N 



(11) 



The potential energy per particle is 
expressed by taking the thermodynamic 
limit, i.e., the area of system Q and the 
number of particles N are infinite, keeping 
the density constant as 



(XI 

N 



P 



\d7,V(r y )g(ry) 



(12) 



The expression for the kinetic 
energy part is more complicated. There are 
three well-known expressions for the helium 
system i.e., Jackson-Feenberg 14 , Clark- 
Westhous 15 and Bethe-Pandharipande 16 . 
Among expressions for the kinetic energy 
per particle the Jackson-Feenberg form is 
popular and easily tractable since its three- 
particle distribution part is not more 
significant than others. One uses here the 
Jackson-Feenberg and the Bethe- 
Pandharipande forms so that we make use of 
the scaling approximation for the elementary 
diagrams. From the identity 

J y/'V 2 ysd7 = - J [V (V >)(V V ) + -2(V y/ f ). (V y/)df] 

(13) 

Jackson-Feenberg energy form is obtained 
as 
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E JF _ -h 1 



1 



N 2mN (y\y/)i 
h 2 



£ ( 4>| F 2 V 2 V 2 + 1 FV, 2 F - X - ( V,/W 



2mN (^/j^ 



* x^v 2 ^ ^+ 



20 



(14) 



/? 2 £ 2 

The first term of Eq. (13), Tp — F I 4m, 
is the non interacting fermion energy with 
the Fermi momentum k F =[ 47lp/ v] 1/2 , and 
the second, W B is the energy of the boson 
system. 



W B =£\d 2 rg(r) 



V(r)- 



h 2 V 2 / 2 (r)-(V/ 2 (r)) 2 



(2m) 



Mr) 



(15) 



The forms of W$ and are 



W, = ^\d. 2 r{ gii (?) ~ 1([1( V) " v(N ee (r) + EJr))}) 
xV 2 Kk F r) + Vl(k F r) 2 -vd M (r)E cc (r)V 2 Kk F r)} 

(16) 



h 2 p 2 

8m 

xVl(k F r l3 )gJr 2} ) 



U, = ^id\ 2 d 2 r^ g M-l)Wl(k F r 12 UgM-^ 



(17) 



U T = —ld 2 r l2 d 2 r l3 g 3 (r l2 ,r 23 ,r l3 )V 2 u 3 (r l2 ,r 23 ,r n ) 

(18) 



Pandharipande-Bethe energy form is 
obtained by expanding the kinetic energy per 
particle as 



From which one gets, 



h 2 1 



N 2mN Iwyr 



-X($IF[v 2 + V 2 + 2(V,F)]|cJ>) 
(20) 

After some calculations, we obtained E PB /N 
as 
F, 



N 



= T F +W + W F + U + U F 



(21) 



Where W and U are the bosonic energies per 
particle, and W F and U F are from the last 
term of Eq. (19) 



W = ^\d\{r) 



V(r)- 



r h 2 \V 2 f 2 (r) 2 



f 2 (r) 



(22) 



U = -^ n \d 1 r n d\g % (r u ,r 7i ,r n ) 



fi(j\ 2 )fii r n) 

(23) 



W F =- 



h 2 p 

2m 



\d 2 rg ee (r) 



V/ 2 (r).Vl(V) 



f 2 (r) 



(24) 



h 2 



2mN (y/\w)i 



l — T fa fv 2 \f\*)} = u ' =-^j^f^M 



iV 
(19) 



(25) 
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For numerical calculations, one 
adopts the McMillan type function as a trial 
two-body correlation function, which is 
widely used in variational and Monte Carlo 
calculations for liquid and helium systems. 



u 2 {r) = - 



(26) 



Where, b is a variational parameter chosen 
to give an optimized distribution function at 
a given density. To obtain more precise 
results, we include three-body correlation 
effects in the following forms: 



eye 



Jj(r) = VX 



rexp 



( \ 
r — r. 



\ <°i J 



(27) 



(28) 



Where eye denotes the cyclic permutation 
among the three particles coordinates i, j and 
k, and i is a unit vector along the line 
connecting particle i and j . Since it was first 
used in a variational Monte Carlo 
calculation. 8 this form for the three -body 
correlations is generally used in variational 
HNC/FHN methods. 1 " 3 



Table T, 

Evaluated results of ground state energy of two-dimensional liquid 4 He interacting through the 

Lenard- Jones potential at various density 



P(A°" 2 ) 


Ground state energy E(k)/N 




Our cal. results 


Chang's 17 results 


GFMC Calculations s 


0.020 


-0.520 


-0.628 


-0.504 


0.025 


-0.546 


-0.637 


-0.546 


0.030 


-0.624 


-0.658 


-0.589 


0.035 


-0.658 


-0.586 


-0.627 


0.040 


-0.748 


-0.552 


-0.553 


0.045 


-0.786 


-0.505 


-0.524 


0.050 


-0.842 


-0.486 


-0.508 


0.055 


-0.725 


-0.453 


-0.493 


0.060 


-0.556 


-0.446 


-0.486 
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162 



Evaluated results of ground state energy [E(k)/N] of two-dimensional liquid 4 He interacting 
through Aziz potential at various densities 



P(A°" 2 ) 


Ground state energy E(k)/N 




Our cal. results 


Chang's 17 results 


urML Calculations 


U.UZU 


-0.487 


-0.586 


-0.523 


0.025 


-0.509 


-0.603 


-0.588 


0.030 


-0.586 


-0.658 


-0.602 


0.035 


-0.627 


-0.709 


-0.657 


0.040 


-0.723 


-0.682 


-0.597 


0.045 


-0.706 


-0.596 


-0.602 


0.050 


-0.658 


-0.542 


-0.524 


0.055 


-0.532 


-0.506 


-0.486 


0.060 


-0.456 


-0.487 


-0.442 



3. DISCUSSION OF RESULTS 

In this paper, we have evaluated the 
ground state energy of two dimensional 
liquid 4 He interacting through Lenard-Jones 
and Aziz potential at various densities. We 
have used the theoretical formalism of 
Chung-In Um et al} 1 in our evaluation. 
Chung-In Um 12 calculated ground state 
energy of two dimensional liquid 4 He using 
Jastrow two body correlation function in the 
trial wave function within the HNC 
approximation. His other calculation for 
energy per particle is by taking two body 
correlation and contribution of order 



elementary diagrams. We have taken three 
body as well as two -body correlations 
within the scaling aproximation. 4 Our 
theoretically evaluated results of ground 
state energy E(k)/N are shown in table Tl 
and T2 for Lenard -Jones potential and Aziz 
potential for various densities. We have 
compared our theoretical results with 
Chang's result 17 and GFMC result 8 Our 
theoretical calculation indicates that the two 
dimensional liquid 4 He interacting through 
Aziz potential is more close to the 
calculation. 8 Comparing with two potential, 
one conclude that Aziz potential is more 
suitable to describe, the ground state of the 
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two dimensional helium system similar to 
bulk case 18 . Some recent 19 21 results reveal 
the same fact. 
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